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Stable counteralignment of a circumbinary disc 
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ABSTRACT 

In general, when gas accretes on to a supermassive black hole binary it is likely to have no 
prior knowledge of the binary angular momentum. Therefore a circumbinary disc forms with 
a random inclination angle 6 to the binary. It is known that for 6 < 90° the disc will coalign 
wrt the binary. If ^ > 90° the disc wholly counteraligns if it satisfies cos < -/d/2/b. where 
/(J and /b are the magnitudes of the disc and binary angular momentum vectors respectively. 
If however > 90° and this criterion is not satisfied the same disc may counteralign its 
inner regions and, on longer timescales, coalign its outer regions. I show that for typical disc 
parameters, describing an accretion event on to a supermassive black hole binary, a misaligned 
circumbinary disc is likely to wholly co- or counter-align with the binary plane. This is 
because the binary angular momentum dominates the disc angular momentum. However with 
extreme parameters (binary mass ratio M2/M1 ^ 1 or binary eccentricity ^ ~ 1) the same disc 
may simultaneously co- and counter-align. It is known that coplanar prograde circumbinary 
discs are stable. I show that coplanar retrograde circumbinary discs are also stable. A chaotic 
accretion event on to an SMBH binary will therefore result in a coplanar circumbinary disc 
that is either prograde or retrograde with respect to the binary plane. 
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1 INTRODUCTION 

Galaxy mergers are commonly thought to be the main mechanism 
driving the coevolution of galaxies and their central supermassive 
black holes (SMBH). An SMBH binary is likely to form in the cen- 
tre of the merged galaxy and subsequently accrete from a circumbi- 
nary disc. It is reasonable to expect that the angular momentum of 
the binary and that of the accreting gas are uncorrelated, and that 
each has no prior knowledge of the other. We can therefore expect 
a random distribution of orientations for such circumbinary discs. 

Understanding the evolution of a misaligned circumbinary 
disc is needed if we are to understand the evolution of SMBH bina- 
ries. In a recent paper (Nixon et al., 2011) we showed that a retro- 
grade circumbinary disc can be very efficient in extracting angular 
momentum from the binary orbit. This may off'er a solution to the 
final parsec problem (Begelman et al. 1980; Milosavljevic & Mer- 
ritt 2001). 

In another recent paper (Nixon, King & Pringle, 2011, here- 
after NKP) we showed that the dominant effect of the binary po- 
tential, on a misaligned circumbinary disc, is to induce radially- 
dependent precessions of the misaligned disc particle orbits. This 
differential precession is known to induce warping of the disc as 
rings of gas dissipate energy through viscosity. The precession van- 
ishes only when the orbit is in the plane of the binary (either pro- 
grade or retrograde) and thus we expect that both prograde and ret- 
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rograde orbits are possible equilibria for the gas. NKP showed that 
the evolution of a misaligned circumbinary disc is formally similar 
to the evolution of a misaligned disc around a spinning compact 
object. In this case the evolution is driven by the Lense-Thirring 
effect (e.g. Bardeen & Petterson 1975; Pringle 1992; Scheuer & 
Feiler 1996; Lodato & Pringle 2006; Nixon & King 2012). This 
implies that the analysis of King et al. (2005), which calculates 
the conditions for co- or counter-alignment, holds for circumbi- 
nary discs as well. The disc co- or counter-aligns depending on 
the magnitudes and directions of 7b and/d, the angular momentum 
of the binary and disc respectively. The whole disc counteraligns 
with the binary if the initial inclination angle of the disc, 6, and the 
magnitudes of the disc and binary angular momentum, and /b 
respectively, satisfy 

cos^<-A. (1) 
z7b 

NKP only considered the zero-frequency (azimuthally sym- 
metric m = 0) term in the binary potential. This is a reasonable 
approach as all other terms induce oscillatory effects which cancel 
out on long timescales. 

It is natural to assume that coaligned discs are stable, how- 
ever in the past counteraligned discs have been incorrectly found 
to be unstable (Scheuer & Feiler, 1996). I therefore use a full three 
dimensional hydrodynamic approach to check the assumptions in 
NKP are valid and to confirm that co- and counter-alignment are 
stable. In section 2 1 discuss the possibility of both co- and counter- 
alignment of the same disc (cf Lodato & Pringle 2006). In section 3 
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I report a simulation of a counteraligning disc and give my conclu- 
sions in section 4. 



2 SIMULTANEOUS CO- AND COUNTER-ALIGNMENT? 

Lodato & Pringle (2006) considered the alignment of a disc and a 
spinning black hole. They showed that for a disc where 6 > njl 
and /d > 2/h, initial counteralignment of the inner disc occurs, 
followed by subsequent coalignment of the outer disc. During the 
alignment of the outer disc the inner disc retains its retrograde na- 
ture and so a warp of significant amplitude is achieved (A^ ~ n). 
This scenario produces a disc which is simultaneously counter- 
aligned (in the inner parts) and coaligned (in the outer parts). For 
the subsequent evolution of such a disc see Nixon, King & Price 
(2012). 

If however the angular momenta are such that < /h then the 
disc can only wholly co- or counter-align and the above process is 
impossible. In this section we discuss whether it is feasible to have 
circumbinary discs such that both co- and counter-alignment are 
simultaneously possible. This possibility is constrained as the mass 
of the disc is limited by self-gravitational collapse and the radius 
of the disc is limited by the gravitational sphere of influence of the 
binary. Therefore there must be a maximum feasible disc angular 
momentum. 

The disc simultaneously co- and counter-aligns if (and only 
if) 6 > nil and the condition (1) does not hold. So for randomly 
aligned discs, simultaneous co- and counter-alignment is only pos- 
sible if /d > 2/b (assuming cos ^ ~ -1). Thus I derive the condition 
for /d ^ 2/b- The disc angular momentum is 



Md ^GMR^ 



(2) 



where Md is the mass in the disc, M is the total binary mass, is a 
characteristic radius for the disc and G is the gravitational constant. 
We note that the definition of in this case is simply the total 
angular momentum in the disc. The exact definition of /d is usually 
more subtle because the disc takes time to communicate its angular 
momentum. 

The angular momentum of a binary with eccentricity e is 



7b = fi-sjGMa(l-e^) 



(3) 



where M is the total mass of the binary and is the reduced mass 
(jjL ^ M2 for M2 <^ Ml). 

So for simultaneous co- and counter-alignment we require 



Md ^JGMRd > 2/1 ^GMa(l-e^). 

Self-gravity limits Md < (H/R)M, so we get 



(H/R)M ^JGMRa > 2/i -yjGMa(l - e^). 



(4) 



(5) 



After a bit of algebra this tells us that we can only get simul- 
taneous CO- and counter-alignment if 



(6) 



For the most optimistic parameters we have a ~ 0. 1 pc and 
the sphere of influence of the binary ~ 10 pc. This gives the LHS 
as < 100. For typical disc thickness (H/R ~ 10~^) it is clear that 
binaries with a low eccentricity (e ^ 0) we cannot get simultane- 
ous CO- and counter-alignment unless the mass ratio is extreme 
(M2/M1 < 10~^). However if the disc is very thick then this is pos- 
sible. For small binary mass-ratios the disc mass may be greater 



than the mass of the secondary, therefore the hydrodynamic drag 
on the binary may become significant (Ivanov et al., 1999). If the 
binary is significantly eccentric, as predicted for binaries that have 
accreted through a retrograde disc (Nixon et al., 2011) then simul- 
taneous CO- and counter-alignment may be possible. I shall return 
to these possibilities in future work. However these arguments sug- 
gest that it is reasonable to assume an SMBH binary must wholly 
CO- or counter-align the disc with the binary plane. 



3 SIMULATION 

To confirm the stable counteralignment of a circumbinary disc I 
perform one simple simulation. I use the SPH code phantom, a low- 
memory, highly efficient SPH code optimised for the study of non- 
self-gravitating problems. This code has performed well in related 
simulations. For example Lodato & Price (2010) simulated warped 
accretion discs and found excellent agreement with the analytical 
work of Ogilvie (1999) on the nature of the internal accretion disc 
torques. 

The implementation of accretion disc or-viscosity (Shakura & 
Sunyaev, 1973) in phantom is described in Lodato & Price (2010). 
Specifically, we use the 'artificial viscosity for a disc' described 
in Sec. 3.2.3 of Lodato & Price (2010), similar to earher SPH ac- 
cretion disc calculations (e.g. Murray, 1996). The main diff'erences 
compared to standard SPH artificial viscosity are that the disc vis- 
cosity is applied to both approaching and receding particles and that 
no switches are used. The implementation used here diff'ers slightly 
from Lodato & Price (2010) in that the term in the signal ve- 
locity is retained in order to prevent particle interpenetration. The 
disc viscosity in phantom was extensively calibrated against the ID 
thin Qf-disc evolution in Lodato & Price (2010) (c.f. Fig. 4 in that 
paper) and the disc scale heights employed here are similar. We use 
Qf^^ = 1 and = 2 which at the employed resolution (see below) 
corresponds to a physical viscosity ass ~ 0-05. Note that the initial 
viscosity is slightly smaller than this, however during the simula- 
tion the disc spreads and particles are accreted and thus the spatial 
resolution decreases which increases the viscosity (see eq. 38 of 
Lodato & Price 2010). 



3.1 Setup 

The simulation has two equal mass sink particles representing the 
binary (of total mass unity in code units), on a circular orbit with 
separation 0.5 (in code units). Initially the circumbinary gas disc is 
flat and composed of 1 million SPH particles, in hydrostatic equi- 
librium from 1.0 to 2.0 in radius, and surface density distribution 
Z oc set up by means of the usual Monte-Carlo technique. 
The vertical hydrostatic equilibrium corresponds to H/R - 0.05 at 
R - \. The equation of state for the gas is isothermal. The disc 
is initially tilted at 170° to the binary plane. Any gas which falls 
within a radius of 0.5 from the binary centre of mass is removed 
as it no longer has any eff'ect of the alignment of the circumbinary 
disc. The disc mass is negligible in comparison to the binary mass 
and thus the gravitational back-reaction of the gas on the binary is 
not included. 



3.2 Geometry 

The binary is taken to orbit in the x-y plane with binary angular 
momentum vector in the z-direction. We define the 'tilt' and the 
'twist' of the disc with respect to the binary using Euler angles (e.g. 
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Figure 1. The tilt angle in the disc as a function of radius. This is plotted 
at six different times corresponding to (in units of the dynamical time at 
R = \)t = {) (black), t = 500 (red), t = 1000 (green), t = 1500 (blue), 
t = 2000 (grey), t = 2500 (purple). Initially the disc is flat with a global tilt 
of 170°. In time the disc spreads and counteraligns wrt the binary plane. 



Figure 2. The twist angle (angle from x-axis to the line of nodes) atR = 1 in 
the disc as a function of time. Initially the line of nodes is in the i/-direction 
and hence the angle is 90° . The binary potential induces precession of the 
gas orbits. The twist is calculated between +180° which generates the 'saw- 
tooth' structure in the plot. 



Bardeen & Petterson 1975; Pringle 1996), where the unit angular 
momentum vector in the disc is described at any radius by 

i = (cos 7 sin sin y sin j3, cos /?) (7) 

with /3 (R, t) the local angle of disc tilt with respect to the z-axis and 
7 {R, t) the local angle of disc twist measured from the x-axis. 

3.3 Stable counteralignment 

As predicted by NKP the dominant effect of the binary on the disc 
is to induce precessions in the gas orbits. As the radial range of the 
disc in this simulation is only a factor of two, the precession rate 
changes little between the inner and outer parts of the disc. How- 
ever there is still a differential precession across the disc, which 
leads to a twist. This causes dissipation between rings of gas and 
so a small amplitude warp in the disc (cf. Fig. 1). The disc angular 
momentum vector precesses around the binary angular momentum 
vector for the duration of the simulation. In Fig. 2 the twist angle 
in the disc (at unit radius) is plotted against time. As the twist angle 
is defined between +180° the plot has a 'saw-tooth' shape, as the 
disc precesses until the twist is -180° which is then equivalent to 
180°. 

Figs. 3 & 4 show the disc structures at various times. Fig. 3 
shows the disc column density viewed face on, and Fig. 4 shows it 
edge on. In Fig. 3 we see the usual viscous spreading of the disc, 
and also the lack of any resonances in the disc. This is expected in a 
counterrotating circumbinary accretion disc (Papaloizou & Pringle 
1977; Nixon et al. 2011). In Fig. 4 we can see the tilt of the disc. 
Initially the disc is tilted by 170°, but then precesses and counter- 
aligns with the binary plane. 

I have also performed simulations (not illustrated) of a mod- 
erately eccentric binary {e = 0.3) and an unequal mass ratio binary 
(M2/M1 = 0.5) which both display the same behaviour. 



4 CONCLUSIONS 

In this paper I have shown that for realistic parameters a circumbi- 
nary disc must usually wholly co- or counter-align with a binary. 



However for extreme mass ratios or high eccentricities the binary 
may be dominated by the disc angular momentum (cf. eq. 6). In 
this case evolution similar to Fig. 9 of Lodato & Pringle (2006) is 
expected with simultaneous co- and counter-alignment of the disc. 

I have shown that a disc with an initial inclination of 170° 
to the binary plane stably counteraligns. Circumbinary discs, with 
Jd <^ Jh and an initial misalignment angle of > 90° counteralign 
with respect to the binary. If the disc angular momentum is not 
negligible. King et al. (2005) showed that the condition for coun- 
teralignment is 

cos^<-A. (8) 

A counteraligned circumbinary disc is efficient at shrinking 
the binary as it directly absorbs negative angular momentum when 
capturing gas into circumprimary or circumsecondary discs (Nixon 
et al., 2011). This interaction increases the binary eccentricity. 
Nixon et al. (2011) show that the timescale to increase the eccen- 
tricity from zero to unity is ~ M2/M where M is the mass inflow 
rate through the retrograde circumbinary disc. Once the eccentric- 
ity is high enough gravitational wave losses will drive the binary to 
coalescence. 
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Figure 4. Edge-on column density rendering of the disc at various times in the simulation. Each component of the binary is represented by a red filled circle. 
Over time the disc spreads and precesses around the binary. The precession induces dissipation in the disc which aligns the disc with the binary plane. 
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